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ABSTRACT 

In recent years there has been considerable activity in using gamma-ray bursts 
as cosmological probes for determining global cosmological parameters comple- 
menting results from type la supernovae and other methods. This requires a 
characteristics of the source to be a standard candle. We show that contrary 
to earlier indications the accumulated data speak against this possibility. An- 
other method would be to use correlation between a distance dependent and a 
distance independent variable to measure distance and determine cosmological 
parameters as is done using Cepheid variables and to some extent Type la su- 
pernovae. Many papers have dealt with the use of so called Amati relation, first 
predicted by Lloyd, Petrosian and Mallozzi, or the Ghirlanda relation for this 
purpose. We have argued that these procedure involve many unjustified assump- 
tions which if not true could invalidate the results. In particular, we point out 
that many evolutionary effects can affect the final outcome. In particular, we 
demonstrate that the existing data from Swift and other earlier satellites show 
that the gamma-ray burst may have undergone luminosity evolution. Similar 
evolution may be present for other variables such as the peak photon energy of 
the total radiated energy. Another out come of our analysis is determination of 
the luminosity function and the comoving rate evolution of gamma-ray bursts 
which does not seem to agree with the cosmic star formation rate. We caution 
however, that the above result are preliminary and includes primarily the effect 
of detection threshold. Other selection effects, perhaps less important than this, 
are also known to be present and must be accounted for. We intend to address 
these issues in future publications. 
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1. INTRODUCTION 

The change in our understanding of gamma-ray bursts (GRBs) in less than a decade 
has been unprecedented. We have gone from groping for ways to determine their distances 
(from solar system to cosmological scales) to attempts to use them as cosmological probes. 
Observations by instruments on board a series of satellites starting with BeppoSAX and 
continuing with HETE, INTEGRAL and Swift, have been the primary source of this change. 
The higher spatial resolution of these instruments has allowed the measurement of redshifts 
of many well-localized GRBs, which has in turn led to several attempts to discover some 
emission characteristics which appears to be a "standard candle" (SC for short), or shows 
a well defined correlation (with a small dispersion) with another distance independent mea- 
surable characteristic. One can use such relations to determine the distances to GRBs in a 
manner analogous to the use of the Cepheid variables. Example of this are the lag- luminosity 
and variability-luminosity relations (Norris ct al 2000, Norris 2002, Fenimorc & Ramirez- 
Ruiz 2000, Reichart et al 2001) which were exploited for determining some cosmological 
aspects of these sources (Lloyd, Fryer & Ramirez-Ruiz 2002, Kocevski & Liang 2006) us- 
ing the methods developed by Efron & Petrosian (1992, 1994, 1999). More recently there 
has been a flurry of activity dealing with the observed relation between the peak energy 
Bp of the vFy spectrum and the total (isotropic) gamma-ray energy output 8\so {Ep oc £^^^, 
7] ~ 0.6) predicted by Lloyd et al. (2000, LPMOO) and established to be the case by Amati et 
al. (2002) (see also Lamb et al. 2004, Attiea et al 2004). Similarly Ghirlanda et al. (2004a) has 
shown a correlation with smaller dispersion between Ep and the beaming corrected energy 
Sj (Ep oc S:i', rj' ~ 0.7) where: 

= £,„ X (^-"'f'-y^ (1) 

and ^jet is the half width of the jet. 

These are followed by many attempt to use these relations for determining cosmological 
parameters (Dai et a/. 2004, Ghirlanda et a/. 2004b, 2005 [Gea05], Friedman & Bloom 2005, 
[FB05]). There are, however, many uncertainties associated with the claimed relations and 
even more with the suggested cosmological tests. The purpose of this paper is to investigate 
the utility of GRBs as cosmological tools either as SCs or via some correlation. In the next 
section we review the past and current status of the first possibility and in §3 we discuss the 
energy-spectrum correlation and whether it can be used for cosmological model parameter 
determination. Finally in §4, we address the question of cosmological luminosity and rate 
density evolution of GRB based on the existing sample with known redshifts. 



2. Standard Candle? 

The simplest method of determining cosmological parameters is through SCs. Type la 
supernovae (SNIa) are a good example of this. But currently their observations are limited to 
relatively nearby universe (redshift z <2). Galaxies and active galactic nuclei (or AGNs), on 
the other hand, can be seen to much higher redshifts (z > 6) but are not good SCs. GRBs 



- 3 - 



are observed to similar redshifts and can be detected to even higher redshifts by current 
instruments. So that if there were SCs they can complement the SNIa results. In general 
GRBs show considerable dispersion in their intrinsic characteristics. The first indication 
that GRBs might be SCs came from Frail et al. (2001) observation showing that for a sample 
of 17 GRBs the dispersion of the distribution of is significantly smaller than that of £^iso- 
The determination of requires a well defined light curve with a distinct steepening. The 
jet opening 6 jet depends primarily on the time of the steepening and the bulk Lorentz factor, 
but its exact value is model dependent and depends also weakly on Siso and the density of 
the background medium 



^,., = 0.101 radian x^^j (-j i^-^ j [w^s) ' 

(2) 

In Figure [1] we show the distribution of £^iso and £^ for 25 pre-Swift GRBs (mostly 
compiled in FB05). As evident, there is little difference between the two distributions (except 
for their mean values) and neither characteristics is anywhere close to being a SC. 

We have calculated the jet angle 9 jet (from equation [2]) for 25 pre- Swift GRBs with 
relatively well defined threak- Assuming a gamma-ray efficiency rj = 0.2 and a value of circum 
burst density estimated from broadband modeling of the lightcurve when available (otherwise 
we use the default value of n = 10 cm~^). 

Unfortunately fewer than expected Swift GRBs have optical light curves and their X- 
ray light curves show considerable structure (several breaks and flaring activity) with several 
GRBs showing no sign of jet-break or beaming (Nousek et al. 2006). This has brought the 
whole idea of jet breaks and calculating £^ into question. The upshot of this is that Swift 
like Siso also has a broad distribution extending over two decades. Thus, any cosmological 
use of GRBs must include the effects of the breadth of the distribution^. 



3. Correlations 

When addressing the correlation between any two variables, one should distinguish 
between a one-to-one relation and a statistical correlation. In general the correlation between 
two variables (say, Ep and £^iso) can be described by a bi-variate distribution ip{Ep, Siso)- If 
this is a separable function, ip{Ep,£iso) = (f){Siso)C{Ep), then the two variables are said to 
be uncorrelated. A correlation is present if some characteristic (say the mean value) of one 
variable depends on the other: e.g. {Ep) = g{Siso)- Only in the absence of dispersion there 



^It should be also noted that there is an observational bias in favor of detecting smaller jet angles 
(i.e. earlier jet-breaks), so that the population as a whole (including those with very late jet-breaks) will 
have even a broader distribution. Note also that a SC Sj means that the Ghirlanda relation would have an 
index rj' 
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Fig. 1. — Distribution of 8\s.o and 8^ for 25 pre-Swift GRBs with evidence for a jet-break 
and beaming. 6"^ distribution is shifted by about two orders of magnitudes compared to 
£\so distribution due to the beaming factor correction. However the dispersion of the two 
distribution are very similar (crg.^^ = 0.68, cr^^ = 0.52) and broad indicating that GRBs 
cannot be assumed to be SCs 
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will be a one-to-one relation; C(-^p) = — g{£iso))- 

In general, the determination of the exact nature of the correlation is complicated by the 
fact that the extant data suffers from many observational selection biases and truncations. 
An obvious bias is that most sample are limited to GRBs with peak fluxes above some 
threshold. There are also biases in the determination of E""^^ (see e.g. Lloyd & Petrosian 
1999; LP99). The methods devised by Efron & Petrosian (1992, 1994) are particularly 
suitable for determination of correlations in such complexly truncated data. 

The first indication of a correlation between the energetics and spectrum of GRBs came 
from Mallozzi et al (1995), who reported a correlation between observed peak flux fp and 
Ep""^. A more comprehensive analysis by LPMOO, using the above mentioned methods, 
showed that a similar correlation also exist between the observed total energy fluence Ftot 
and E!^""^. Both these quantities depend on the redshift z = Z — 1; 

Ftot = SiJ{A7TdlZ), and E;''' = Ep/Z, (3) 
Here dm is the metric distance, and for a flat universe 

dm{Z) = {c/Ho) dZ\n{Z'))-^l\ with fi(Z) = p(Z)/po, (4) 

describing the evolution of the total energy density p{z) of all substance (visible and dark 
matter, radiation, dark energy or the cosmological constant). LPMOO also showed that the 
correlation expected from these interrelationships is not sufficient to account for the observed 
correlation, and that there must be an intrinsic correlation between Ep and £^iso- Without 
knowledge of redshifts LPMOO predicted the relation Eiso oc E^'^ which is very close 0.5, 
which is very similar to the so-called Amati relation obtained for GRBs with known redshifts. 
However, it should be emphasized that the LPMOO result implies a statistical correlation and 
not a one-to-one relation needed for using GRBs as a reliable distance indicator. Nakar & 
Piran (2004, 2005) and Band & Preece (2005) have shown convincingly that the claimed 
tight one-to-one relations cannot be valid for all GRBs. We believe that the small dispersion 
seen in GRBs with known redshifts is due to selection effects arising in the localization and 
redshift determination processes: e.g., these GRBs may represent the upper envelope of the 
distribution. A recent analysis by Ghirlanda et al. (2005) using pseudo-redshift shows a much 
broader dispersion (as in LPMOO). The claimed tighter Ghirlanda relation, could be due to 
additional correlation between the jet opening angle 6'jet and £^iso, Ep, or both. However, as 
mentioned above the picture of jet break and measurements of 9jet and £^ is a confusing 
state in view of Swift observation. 

We have reanalyzed the existing data and determined the parameters of the Amati and 
Ghirlanda relations. In Figure [2] we show the the £^iso (and £^ excluding some outliers) vs Ep 
for all GRBs with known redshifts (and 6'jet)- We compute best power law fit for both these 
correlations and we describe the dispersion around it by the standard deviation. For the 
Ep — £iso we find: Ep oc E^^^, rj ~ 0.328 ± 0.036 and aiso = 0.286. For Ep — £^ correlation we 
find: Ep (X SJ^ , rj' ^ 0.555 ± 0.089 and = 0.209. We find that additional data has reduced 
the significance of the correlations or has increased the dispersions (compare our values of 
aiso ~ .... and ~ .... obtained by Amati et a/, and Ghirlanda et al. ). This is contrary to 
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what one would expect for a sample with a true correlation. From this we conclude that, as 
predicted by LPMOO, there is a strong correlation between Ep and £^iso (or S^), but for the 
population GRBs as a whole both variables have a broad distribution and most GRBs do 
not obey the tight relations claimed earlier. 



3.1. Correlations and Cosmology 

Attempts to use observations of extragalactic sources for cosmological studies have 
shown us that extreme care is required. All observational biases must be accounted for 
and theoretical ideas tested self-consistently, avoiding circular arguments. This is especially 
true for GRBs at this stage of our ignorance about the basic processes involved in their 
creation, energizing, particle acceleration and radiation production. Here we outline some 
of the difficulties and how one may address and possibly overcome them. 

Let us assume that there exists a one-to-one but unknown relation between £^iso and Ep, 
^iso = ^ofiEp/Eo), and that we have a measure of Ftot and z. Here £q and Eq are some 
constants, and for convenience we have defined f{x) which is the inverse of the function g 
introduced above. The 



From equations ([3]) and (jlj) we can write 



For general equations of state P = Wip, Q{Z) = The aim of any cosmological 

test is to determine the values of different Qi and their evolutions {e.g. changes in Wi) . If we 
make the somewhat questionable assumption of complete absence of cosmological evolutions 
of Siso,Ep and the function f{x), then this equation involves two unknown functions Q{z) 
and f{x). In principle, if the forms of these functions are known, then one can rely on some 
kind of minimum method to determine the parameters of both functions, assuming that 
there is sufficient data to overcome the degeneracies inherent in dealing with large number of 
parameters. By now the parametrization of Q{Z) has become standard. However, the form 
and parameter values of f{x) is based on poorly understood data and theory, and currently 
requires an assumed cosmological model. Using the form {e.g. the power law used by Gea05) 
derived based on an assumed cosmological model to carry out such a test is strictly speaking 
circular. (It is even more circular to fix the value of parameters, in this case the index rj, 
obtained in one cosmological model to test others as done by Dai et al. 2004). Even though 
different models yield results with small differences, this does not justify the use of circular 
logic. The differences sought in the final test using equation [5] will be of the same order. The 
situation is even more difficult because as stressed above the correlation is not a simple one- 
to-one relation but is a statistical one. Finally, the most important unknown which plagues 
all cosmological tests using discrete sources is the possibility of the existence of an a priory 
unknown evolution in one or all of the relevant characteristics. For example, the intrinsic 
luminosity Lj^o might suffer large evolution which we refer to as luminosity evolution. The 
value of Ep can also be subject to selection effects, or the correlation function f(x) may evolve 
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Siso — Epl 8^ — Ep correlations 
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Fig. 2. — Ep — Siso and Ep — correlations. The 43 gray circles are all the bursts from 
our sample that had good enough spectral observations to find the energy peak of the uF^ 
spectrum, and the 21 black diamonds are a subset of those bursts with a jet break found in 
their optical lightcurve. Sohd fines are the best fit we find for the two correlations {Ep oc S^^^, 
77 ~ 0.328 ± 0.036 and Ep oc rj' ~ 0.555 ± 0.089) and dashed fines are the best fit found 
by Amati et al. 2003 with 20 bursts {Ep oc E^^f) and Ghirlanda et al. 2004 with 16 bursts 
(^pocffo). 
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with redshift, i.e. t] = r]{z). For such a general case we are dealing with 4 unknown ??? of 
the two above. Moreover the rate function of GRBs most likely is not a constant and can 
influence the result s with a broad distribution. We address some of these questions now. 



4. GRB Evolutions 

For a better understanding of GRBs themselves and the possibility of their use for cos- 
mological tests we need to know whether characteristics such as £iso, 6'jct, Ep, the correlation 
function f{x) and the occurrence rate pgrb (number of GRBs per unit co- moving volume 
and time) change with time or Z. For example to use the Ep — £^iso correlation for cosmolog- 
ical purposes, one need to first establish the existence of the correlation and determine its 
form locally (low redshift). One then has to rely on a theory or non-circular observations to 
show that either this relation does not evolve or if it does how it evolves. The existing GRB 
data is not sufficient for such a test. In factthere seems to be some evidence that there is 
evolution. Lie??? has shown by subdividing the data into 4 z-bins, he obtained different 
index r) which change significantly, rendering previous use of this relation for cosmological 
test invalid. This emphasizes the need for a solid understanding of the evolution of all GRB 
characteristics. Two of the most important characteristics are the energy generation £iso and 
the rate of GRBs. These are also two characteristics which can be determined more readily 
and ??? with higher uncertainty. In what follows we address these two questions. We will 
use all GRBs with known Z irrespective of whether we know the jet angle ^jct because this 
gives us a larger sample and because in view of new Swift observations (Nousek et al. 2006) 
the determination of the latter does not seem to be straightforward. Also since it is often 
easier to determine the peak flux fp rather than the fluence threshold, in what follows we 
will use the peak bolometric luminosity Lp = Aird'^Z'^fp instead of Siso = J L{t)dt. 



4.1. Evolution with Pseudoredshifts 

Before considering GRBs with known redshifts we briefly mention that there has been 
two indications of strong evolutionary trends from use of pseudo redshifts based on the so- 
called luminosity- variability and lag-luminosity correlations (Lloyd et al. 2002, Kocevski & 
Liang 2006) using the methods developed by Efron & Petrosian (1992, 1994). These works 
show existence of a relatively strong luminosity evolution L{z) = LqZ"' {a = 1.4 ±0.5, 1.7 ± 
0.3) from which one can determine a GRB formation rate which also varies with redshifts 
and can be compared with other cosmological rates such as the star formation rate. 



- 9 - 



4.2. Evolution with Measured Redshifts 

4.2.1. Description of the Data 

We have compiled the most complete list of GRBs with known redshift. Since the 
launch of the Swift satellite, this list has become significantly larger. We include only bursts 
with good redshift determination meaning that GRBs with only upper or lower limits on 
their redshift are not in our sample. On total, our sample contains 86 bursts, triggered 
by 4 different instruments: BATSE on board CGRO (7 bursts), BeppoSAX (14), HETE- 
2 (13), and Swift (52). For each burst we collected fiuence and peak flux in the energy 
bandpass of the triggering instrument, as well as the duration of the burst. When available 
we have also collected spectral information namely the parameters that define the Band 
function; the energy peak {E°^'^) as well as the low [a) and high (/3) energy indexes of the 
vFy spectrum. When a good spectral analysis was not available, we took as default values 
the mean of the BATSE distributions based on large number of bursts: < a >= —1.0, 
< f3 >= —2.3 and < E"^"^ >= 250 keV. For non-Swift bursts, all this information was mostly 
extracted from FB05 data set to which we added results from recent spectral analysis released 
in GCN. For Swift GRBs, red shift, duration, fiuence and peak flux were co mpiled from 
the Swift Information webpage (jhttp: //swift. gsfc.nasa.gov/docs/swift/archive|) and spectral 
information have been retrieved from GCN releases and we have also looked at spectral 
analysis ourselves for some of them. We assumed the following cosmological model: 



= 0.3, = 0.7, Ho = 70 km s^^ Mpc^^ 

in order to determine the intrinsic properties {e.g. £iso see eq. [3]). £iso is here calculated for 
a rest-frame bandpass [20,2000] keV. Note that K-correction due to the shift of the photons 
into the instrument bandpass has been properly taken into account for the £iso calculation 
(see Bloom et a/. 2001). From this, we calculate the average isotropic-equivalent luminosity 
as: 



90 



when T90 is the duration of the burst that includes 90% of the total counts. 

Because different instruments have been used to collect this information, the sample 
is very heterogeneous and suffers from various selection and truncation effects that vary 
from burst to burst. The most simple of these truncation effects is due to the limiting 
sensitivity of the instruments. A GRB trigger will occur when the peak flux of the burst 
exceeds the average background variation by a few sigmas (depending on the setting of the 
instrument). In an attempt to carefully take into account this effect into our study, we used 
the analysis carried out by Lamb et al. 2005 for pre- S'w/i instruments. In this analysis, they 
computed the sensitivity for each instrument depending on the spectral parameters of the 
bursts. Therefore, for each specific burst of our data set, from its spectral parameters it is 
possible to determine the limiting photon flux of the instrument (for specific GRB with its 
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Fig. 3. — Isotropic average luminosity versus redshift for all bursts in our sample (86). Dif- 
ferent symbols represent bursts observed by different instruments: BATSE(7), Beppo-SAX 
(14), HETE-2 (13), Swift (52). For all non- Swift hurst, a vertical line is plotted representing 
the range of isotropic luminosity in which the burst would still have been observable by the 
instrument keeping all its others parameters fixed. Using the work of Lamb & al. 2005, the 
limiting luminosity is taken to be only dependent on the energy peak Ep of the bursts. For 
Swift bursts, a conservative threshold flux of 0.8 ergs s~^ cm~^ has been chosen. This limit 
is shown as a dashed line. 
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specific spectral parameters). From these, we can easily compute the limiting peak flux fp^im 
for our burst (assuming the Band function for our spectrum). Finally, we can determine the 
detection threshold of the observed energy fluence Fobs- This lower limit Fobs,um. is obtained 
via the simple relationship (Lee & Petrosian 1996): 

Fobs fp 
Fobs,lim fp,lim 

Using the same reasoning we can obtain the limiting values for the intrinsic quantities mean- 
ing the intrinsic values that a given burst needs to have in order to be detected: 

^iso LigQ fp 

^iso,lim LigQiiYn fp,lim 

Those limiting average luminosities for each bursts of our sample are represented in Figure 
[3l This analysis was not carried out for BAT instrument on board Swift therefore we used 
a conservative threshold of 0.8 erg s~^ cm~^ for all of Swift bursts. 



4-2.2. Analysis and results 

We now describe our determination of luminosity and density rate evolution of the 
parent population of our GRB sample. Our analysis is based on the work done by Efron & 
Petrosian (1992, 1994). We refer the reader to these two papers for details. We will here 
simply describe the most important steps of the analysis and what it allows us to infer on 
our data sample. This method has been developed in order to take into account effects of 
data truncation and selection bias on a heterogeneous sample from different instruments with 
different sensitivities as described above and shown in Figure [31 The method corrects for 
this bias by applying a proper rankings to different subset of our sample. The flrst step is to 
compute the degree of correlation between the isotropic luminosity and redshift. For that we 
use the specialized version of Kendell's r statistics. The parameter r represents the degree 
of correlation found for the entire sample with proper accounting for the data truncation, 
r = means no correlation is found between the two parameters being inspected (luminosity 
and redshift in our case). Any other speciflc value Tq implies presence of a correlation with a 
signiflcance of tqct. With this statistic method in place, we can calculate the parametrization 
that best describe the luminosity evolution. To establish a functional form of the luminosity- 
redshift correlation, we assume a power law luminosity evolution: L{z) = LoZ'^. We then 
remove this dependency from the observed luminosity: L' L observed/ {1 + z)" and calculate 
the Kendell's r statistics as a function of a. Figure H] shows the variation of r with a. 

Once the parametric form for the luminosity evolution have been determined, this non- 
parametric maximum likelihood techniques can be used to determine the cumulative distri- 
bution for luminosity and redshift (see Efron & Petrosian 1994) say ^{L) and cr(2;), which 
gives the relative number of bursts under a certain redshift z. From this last function, we can 



easily draw the comoving rate density n{z), which is the number of GRB per unit comoving 
volume and unit time: 



n{Z) 



da{Z) Z 
dZ dV/dZ 



(6) 



where the factor Z is to take the time dilatation into account. Note that this method 
do not provide any constrain on the normalization of any of the quantity mentioned above. 
Normalization will therefore be set arbitrarily on all our figures representing these functions. 

We find a 3.68 a evidence for luminosity evolution (see the r value at the onset of Figure 
Hlwhen a = 0). From this figure, we can also infer that a = 2.21 for value obtained when 
r = gives the best description of the luminosity evolution for the assumed form has a one 
sigma range of [1.75,2.74]. Constrain on the a parameter is not very tight off course since 
the size of our sample is still limited. While current satellites accumulate more data, we will 
be able to increase our data set and further constrain this parameter in the future. 
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Fig. 4. — Variation of the r parameter with the power law index a of the luminosity evolution, 
r = means no correlation which gives the best value of a = 2.21 for the assumed power 
law form with a one sigma range of 1.75 to 2.74. 
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Fig. 5. — The cumulative luminosity distribution $(!/) (left panel) and the cumulative 
redshift distribution a{z) (right panel). Fit to the cumulative density rate is represented 
with a point line. 

The cumulative functions are both shown in Figure [5] and the estimated comoving 
density rate is shown by the jagged curve in Figure 6. Most of the high frequency variation 
is not real and is due to taking the derivative of a noisy curve {a{Z)). The dashed line 
was obtained by fitting the cumulative density distribution by the following parametrized 
function: 

with the following values for the parameters: Zq = 1.8, pi = 7.1, and p2 = 0.95 

These results are still very preliminary as more data become available, accuracy of the 
density function will increase constraining further the evolution rate of long bursts. By 
tackling this problem for the first time we hope to set the ground for further analysis in the 
future. 

The behavior of the comoving density rate for our sample of long bursts is quite peculiar 
with a significant rate increase happening at low redshifts. This effect might be due to some 
selection effects that we have not included in our analysis. For instance, it might be a 
consequence of the fact that instruments detect more easily low-redshifts host galaxies and 
therefore create a bias toward low redshifts GRBs. Another interesting feature is the steady 
increase we obtain in the GRB rate at high redshifts {z > 3). Figure 7 compares the 
estimated comoving rate evolution with different models of Star Formation Rates (SFRs I, 
II, III). 

For comparison with Star Formation history, we used three different models taken from 
the literature: 
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Fig. 6. — The comoving rate density n{z). The dashed line was obtained by fitting the 
cumulative density distribution by the parametrized smooth function of equation [71 We also 
show comparison of the density rate (from Figure [6]) with three different SFR scenarios taken 
from literature. No SFR scenario seems to match the density rate deduced from our analysis. 
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Fig. 7. — Comparison of the comoving density rate evolution of the total sample with that 
of several sub-sample where we impose three different luminosity thresholds: Lj^o > 10''^, 
> 10^°, and > 8 X 10^° ergs.s"^. We also looked at a low luminosity population where we 
imposed of maximal luminosity of 8 x 10^° ergs.s"^. Each sub-sample is subject to the same 
analysis and has provided different luminosity evolution as evident from the different values 
of a. As expected, the rate at low redshift decreases with increasing values of threshold. 
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- Steidel et al. 1999: 



h = 0.16h 



-3 




■70 



e3-4^ + 22 



- Porciani & Madau 2000: 



n = 0.22h' 



Mgyr-^Mpc 
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- Cole et a/. 2001: 



(a + bz)h7o 



M0yr"^Mpc 



-3 



(10) 



n = 



1 + {z/cY 



with {a,b,c,d) = (0.0166,0.1848,1.9474,2.6316) 

As evident, no SFR scenario seems to match the density rate evolution deduced from 
our analysis, specially at low redshift. How much of this difference is real and how much is 
due to other selection effects that we have not quantified is unclear. Because of increasing 
difficulty of identifying the host galaxy with increasing redshift one would expect some bias 
against detection of high redshift bursts. But the largest densities for SFR seems to be in 
the intermediate redshift range. 

An other possibility is that there may exist subclasses of GRBs such as low or high 
luminosity classes. In order to test this eventuality we have defined several subsets of our 
total GRB sample carried out the above analysis for each subsamples, determining a new 
luminosity evolution (a new a) and then proceeding to obtain n{z) from the smooth function 
fitting a{z). We impose different luminosity thresholds for the different subsamples. Three 
different threshold have been chosen: Lj^o > 10"^^ ergs.s"^, Liso > 10^° ergs.s"^, and Liso > 
8x 10^° ergs.s"^. We also looked at a low luminosity population where we imposed of maximal 
luminosity of 10^° ergs.s"^. Figures [6] and [7] compare the new rates with that of the total 
sample. As expected high luminosity samples contribute less to the rate at low redshifts. 
But the general trend and the differences with SFR are essentially still present. However the 
method and framework we presented would be a very valuable tool when enough data has 
been accumulated. 



We have considered GRBs as cosmological tools. We find that GRBs are not SCs 
and the correlations found so far are statistical in nature and too broad to be very useful for 
cosmological model parameter determination. In addition we have shown that there is strong 
evidence for evolution of the peak luminosity of GRBs. This indicate a very likely possibility 
that £^iso and S^/ may also have undergone comparable changes. There may also be evolution 
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of Ep, 6jct or other relevant characteristics. This makes use of GRBs as cosmological tools 
more difficult. 

We may therefore ask is this process hopeless? Strictly speaking the answer is no. Some 
broad brush conclusion can already be reached. For example, one can test the relative merits 
of different forms for f{x). As shown by FB05 the SC assumption (0(£^iso) — (^(^iso — ^o)) 
gives unacceptable fit to essentially all cosmological models, but the use of the power-law 
form agrees with Lemaitre type model (FB05, Gea05) with a relatively long quasi-static 
phase (refereed to as a loitering model). Such models, which were in vogue some time ago 
(see Petrosian 1974), are currently unacceptable because of their low (baryonic plus dark) 
matter density and large curvature. This indicates that the form of the correlation and/or 
other assumptions {e.g. no evolution) are not correct. In this paper, we have shown that 
there is strong evolution of luminosity and £^iso- Lei et al. have also shown that the form of 
the Ep — £iso correlation may evolve. These are tentative results and more data are required 
to determine these evolution trend and their meaning. 

The relevant variables in addition to redshift are E°^' = Ep/{l + z) and Ftot, determi- 
nation of which requires a good description of the total spectrum. We need to know the 
observational selection biases for all the variables and their parameters, and use accurate 
statistical methods to account for the biases and data truncations. From these one can learn 
about the distributions of S^so and Ep and their correlations. 

Obviously Swift observations will be extremely helpful and eventually may provide the 
data required for this complex task. In the near future, however, from analysis of the in- 
coming and archived data we will (and need to) first learn more about the nature of the 
GRBs than cosmological models. Eventually we may have enough information to construct 
a well defined "SC" , which can be used for global cosmological tests as is done using type 
la supernovae. The immediate situation may be more analogous to galaxies where the cos- 
mological tests are rendered complicated because of the multivariate situation and broad 
distributions of the relevant variables. Consequently, over the years the focus of activity has 
shifted from the determination of the few global cosmological parameters to the investigation 
of structure formation, the building process of the black holes and the star formation rate 
(SFR). Similarly, we expect that from investigation of GRBs we will learn about the evolu- 
tion, distributions and correlations of their intrinsic characteristics, and the relationship of 
these with the evolutionary rates of other cosmological sources and the formation rates of 
stars, supernovae and black holes. 

In summary, on the long run cosmological test with GRBs may be possible, either carried 
out with clever statistical methods, or by identification of a subclass of "SCs" . On a shorter 
time scale, we need to learn more about the intrinsic characteristics of GRBs, and provide 
a reasonable theoretical interpretation for them and their cosmological evolution. As an 
example, as shown here, one outcome of our analysis is the determination of the evolution 
rate of GRBs. We have shown that for all GRBs with known redshifts this rate appears to be 
different form the SFR. These differences seems to be present with different subdivision of the 
sample and may be consequences of other selection biases not included in our analysis. But 
the difference is not what one would expect from some possible observational selection...?? 
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